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22 ■ ABSTRACT 

>• . We have undertaken an X-ray survey of a far-infrared flux limited sample of seven 

^ I nearby edge-on starburst galaxies. The data are presented in paper 1 of this series 

^ I by Dahlem, Weaver, & Heckman. Here, we examine the two X-ray-brightest sample 

members NGC 253 and M 82 in a self-consistent manner, taking account of the spatial 
distribution of the X-ray emission in choosing our spectral models. Both galaxies require 
at least three components to model the 0.1 to 10.0 keV spectrum, but the modeling 
is by no means unique. There is significant X-ray absorption in the disk of NGC 253. 
When this is accounted for in spectral fits to the emission between ~ 0.2 and 2.0 keV, 
we find that multi-temperature thermal plasma models with significant underlying soft 
X-ray absorption are more consistent with the imaging data than single-temperature 
models with highly subsolar abundances or models with minimal absorption and non- 
equilibrium thermal ionization conditions, as have been proposed by others. Our three- 
component models do not require absolute abundances that are inconsistent with solar 
values or unusually supersolar ratios of the a-burning elements with respect to Fe (as 
claimed previously). We conclude that with current data, the technique of measuring 
abundances in starburst galaxies via X-ray spectral modeling is highly uncertain. Much 
improved spatial/spectral resolution and increased sensitivity are required. 
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Based on the point-like nature of much of the X-ray emission in the PSPC hard- 
band image of NGC 253, we suggest that a significant fraction of the "extended" X-ray 
emission in the 3 — 10 keV band seen along the disk of the galaxy with ASCA and 
BeppoSAX (Cappi et al.) is comprised of discrete sources in the disk, as opposed to 
purely diffuse, hot gas. If a significant amount of the hard X-ray emission were due to 
unresolved point sources with weak Fe Ka emission (e.g.. X-ray binaries), this could 
explain the low Fe abundances of ~ 1/4 solar derived by Cappi et al. for pure thermal 
models. 

Subject headings: galaxies: individual (NGC 253; NGC 3034 = M82) — galaxies: star- 
burst — galaxies: intergalactic medium — galaxies: evolution — X-rays: galaxies 
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1. Introduction 

NGC 253 and M 82 are the best-studied star- 
burst galaxies in the X-ray band. Their intri- 
cate nature lias been revealed in observations 
with Rosat and Einstein. M82 has a complex, 
central star-forming region that contains unre- 
solved point sources, diffuse disk emission, and 
an X-ray "cone" pointed along the minor axis of 
the galaxy that is consistent with the outflow of 
gas in a partly confined jet (Bregman, Schulman, 
& Tomisaka 1995; Watson, Stanger k, Griffiths 
1984). NGC 253 has distinct components within 
its central region (Fabbiano &l Trinchieri 1984) 
in the form of point sources in the disk, diffuse 
disk emission, unresolved nuclear emission, and 
an extension to the south of the nucleus associ- 
ated with an outflow, visible e.g., in Ha emis- 
sion (Watson, Stanger, k Griffiths 1984). Both 
M 82 and NGC 253 possess X-ray emission that 
extends 6' to 9' outside the plane of the galaxy 
(Schaaf et al. 1989; Fabbiano 1982) with prop- 
erties consistent with a starburst-driven outflow 
(Dahlem, Weaver, and Heckman 1998; hereafter 
paper 1). 

The marked complexity of these objects has 
resulted in contradictory results from instruments 
that cover different bandpasses and different spa- 
tial regions. For example, while the hard X-ray 
detectors onboard EXOSAT and Gin^a measured 
thermal temperatures of kT ~ 6 to 9 keV (Ohashi 
et al. 1990, Schaaf et al. 1989), soft X-ray de- 
tectors such as Einstein, measured significantly 
lower temperatures (Fabbiano 1982). Recent ex- 
periments with a more complete energy coverage 
have at least partly resolved this duality. Obser- 
vations with the Broad Band X-ray Telescope, 
the Advanced Satellite for Cosmology and As- 
trophysics {ASCA), and BeppoSAX show that 
multiple spectral components are necessary to 
model the data between 0.5 and 10 kcV (Pctre 
1993, Ptak et al. 1997 [hereafter P97], Moran & 
Lehnert 1996 [hereafter ML96], Tsuru et al. 1997 
[hereafter T97]). What is not yet clear is whether 



all of these components are required to be ther- 
mal in nature (Ohashi et al. 1990, P97, Cappi et 
al. 1999 [hereafter C99]). 

A clear understanding of the X-ray properties 
of starburst galaxies requires observations with 
a large energy coverage and good spatial resolu- 
tion. The best such dataset currently available 
is contained in the Rosat and ASCA archives. 
M 82 and NGC 253 have been the subject of at 
least six recent papers that focus on the archival 
data (ML96, P97, T97, Strickland et al. 1997, 
Vogler and Pietsch 1998 [hereafter VP98], paper 
1). Results at higher energies are also becoming 
known from BeppoSAX (C99). Based on the in- 
nate complexity of M82 and NGC 253, it is not 
surprising that the interpretations of these data 
by various authors differ. ML96 find a complex 
absorption structure in the core of M 82 and spec- 
tral data consistent with solar abundance ratios. 
On the other hand, P97 and T97 find subsolar 
abundances as well as an overabundance of the 
a-burning elements (Mg, Si, S, Ne) with respect 
to Fe, which suggests that Type II supernova 
remnants dominate the chemical enrichment of 
the X-ray emitting plasma. Since superwinds 
provide a method for enriching the interstellar 
medium (ISM), these results have far-reaching 
implications for studies of supernova rates and 
outflow models of the enrichment of the ISM. It 
is thus critical to establish how reliable the cur- 
rent results are. 

2. The Scope of This Investigation 

The purpose of this paper, which is part of a 
larger study of a sample of seven edge-on star- 
burst galaxies (paper 1), is to thoroughly exam- 
ine the available X-ray data of NGC 253 and M 82 
with the goal of finding a self- consistent spatial 
and spectral model description of the observed X- 
ray emission between 0.1 and 10 keV. Combin- 
ing the ASCA and Rosat datasets is a powerful 
spectroscopy tool, because all emission compo- 
nents can be modeled at one time. The Rosat 
spectra, after excision of point sources, are best- 
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suited to determine the spectral properties of the 
~ 10^ — 10^ K gas emitting the softest X-rays, 
while the hard X-ray (> 1 — 2 keV) components 
inferred from these spectra are subsequently con- 
strained based on the ASCA fits. This pa- 
per addresses a very fundamental physical prob- 
lem, namely the determination of the appropri- 
ate composite spectral model for starburst galax- 
ies, despite, as we will show, an ambiguity in the 
minimum space. Once the most likely physi- 
cal models are determined, similar models can be 
applied to starburst galaxies with poorer quality 
data (for example, fainter galaxies or those at 
high redshift). 

3. Spectral Analysis 

3.1. A Note on the Spectral Modeling 

There is strong evidence that most of the 
soft X-ray emission in NGC 253 and M 82 arises 
from thermal processes. To model the soft X- 
ray emission (< 2 keV) we consider Raymond- 
Smith plasma models (RS; Raymond & Smith 
1991) and Mewe-Kaastra plasma models with Fe 
L calculations by Liedahl (MEKAL; Mewe, Gro- 
nenschild &: van den Oord [1985], Mewe, Lemen 
k van den Oord [1986], Kaastra [1992]). Typical 
thermal temperatures might be similar to those 
found for diffuse gas in normal galaxies, which 
fall near a few tenths of a keV. The X-ray spec- 
trum from such a plasma is rich in line emission 
near 1 keV (Figure 1). Current data are incon- 
clusive as to whether the hard X-ray emission is 
of thermal origin (Ohashi et al. 1990, P97, T97, 
but see C99). The ASCA data are not particu- 
larly sensitive to the Fe K line emission expected 
for a thermal plasma (P97) , and so for simplicity, 
we adopt bremsstrahlung or power law models to 
describe the ^■2-10 keV continuum. 

Our results are tabulated in Tables 1-8. Each 
fit is assigned one or a combination of the follow- 
ing codes for the model used, together with the 
number of the fit if more than one fit is listed. 
The models consist of: 



• "P" - a power law, 

• "B" - a thermal Bremsstrahlung spectrum, 

• "M" - a MEKAL plasma model, 

• "R" - a RS plasma model. 

All models include absorption due to our Gal- 
axy and additional absorption where necessary. 

3.2. The Definition of Different Emission 
Components 

Our imaging analysis and techniques for spec- 
tral extraction are discussed in paper 1. The 

emission components within each galaxy are spa- 
tially identified in the Rosat HRI and PSPC im- 
ages as follows: 

1. NGC 253: 

At a limiting flux level of ~ 6 x 10~^^ 
ergs cm~^ s~^, we identify eighteen point 
sources in the HRI image (Figure 2) within 
a ~ 12' radius of the galaxy core. Based on 
the smoothed and point-source-subtracted 
PSPC images (paper 1), we divide the gal- 
axy into "halo" , "disk" , and "core" regions 
(Figure 3). To examine the band of ab- 
sorption apparent in the 0.25 keV PSPC 
image (paper 1, Figure 12), we also divide 
the elliptical region that delineates the disk 
equally into northern and southern regions 
(each 1/2 of the NE-SW ellipse in Figure 3). 
Compact sources are excised and classified 
as "hard" or "soft" by comparing counts in 
the 0.25, 0.75, and 1.5 keV bands. Their 
spectra are then co-added^ (§ 4.1). 

ASCA observed NGC 253 with the SIS op- 
erating in 4-CCD mode, which provides a 
22' X 22' field of view. NGC 253 is approxi- 
mately 25' X 12' in size and so the galaxy is 

^The hard compact sources are NGC253:DWH 1, 2, 4, 6, 
9, 12, 13, 14, and 15 in the PSPC image; the soft sources 
are NGC253:DWH 3, 16, and 17 (paper 1). 
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contained almost entirely on two of the four 
SIS chips, with only about 15% of the halo 
flux lost off the edges of the chips. For our 
joint PSPC and ASCA analysis, spectra are 
extracted from 10'5 x 15'0 rectangular re- 
gions in the SIS and 15' diameter circular 
regions in the GIS. This technique includes 
as much of the halo flux as possible. For 
a direct comparison of our ASCA results 
with those of P97 (§ 4.2), wc alternatively 
use circular regions identical to theirs (10' 
in diameter). SIS background is accumu- 
lated from the non-source chips. 

2. M82: 

Based on the Rosat images, we divide the 
galaxy into halo and core regions. Of the 
17 point-like sources identified in the PSPC 
image (paper 1, Figure 13), none are lo- 
cated within the galaxy. These sources are 
sorted by hardness ratio and co-added as 
above for spectral fits^. 

4. Spectroscopic Results for NGC 253 
4.1. Modeling the PSPC Spectra 

From our PSPC imaging analysis (paper 1), 
we label the predominant spatial/spectral sources 
as (1) the diffuse gaseous halo, which we identify 
from the images with emission at low X-ray ener- 
gies (0.25 keV band), (2) the galaxy disk, which 
we identify from the images primarily with emis- 
sion at medium X-ray energies (0.75 keV band), 
and point sources (including the core), which we 
identify from the images with emission at high 
X-ray energies (1.5 keV band). 

1. The Diffuse Halo Emission: 

To describe the halo emission we examine 

one- and two-temperature thermal models. 
The simplest case is a single-temperature 

^The soft compact sources are M 82:DWH 1, 2, 8, 9, 10 and 
15; all others have hard spectra (paper 1). 



MEKAL plasma model with Galactic ab- 
sorption (A^H = 0.9x10^*^ cm^^) and so- 
lar abundances. Such a model provides a 
very poor fit to the data (x^ = 96.1 for 
28 d.o.f., model Ml, Table 1). When the 
abundance is a free parameter, the fit im- 
proves dramatically (x^ = 20.1 for 27 d.o.f., 
model M2), but the implied abundance of 
Z = 0.01 Zr:; is extremely low. We cannot 
be sure that this model provides a correct 
physical description of the data. Indeed, 
when we consider the ASCA data, our 
fits suggest generally higher abundances for 
thermal models. Given the limited band- 
pass and poor energy resolution of the 
PSPC, it is possible that continuum flux 
from a warmer gas can contribute to mimic 
abundances that are less than predicted in 
the cooler gas. 

When the data are fitted with a two-component 
thermal plasma model, the abundances can 
be as large as solar but with a fairly large 
uncertainty. If we assume an abundance 
of Z = 0.2 Zq (Table 1, model MM) then 
kTi = 0.141°:°^ keV and kTa = 0.61±0.36 
keV (x^=19.6 for 26 d.o.f.). The observed 
and absorption-corrected 0.1—2.0 keV fluxes 
for this case are 1.24 and 1.89 x 10~^^ ergs 
cm~^ s^^, respectively, with relative contri- 
butions of 40% and 60% from the cool and 
warm components. This ratio is consistent 
with the ratio of extended halo emission 
in the 0.25 keV and 0.75 keV PSPC maps 
(paper 1), and so we are confident that our 
two-component model is appropriate even 
though there is a large uncertainty on the 
abundance. 

2. The Diffuse Disk Emission: 

For the disk emission, a power-law (model 
P) or a MEKAL plasma (model M) provide 
reasonable fits.^ The large photon index of 

^When dealing with the inner regions of the galaxy, the 
choice of background is important. Poorer fits result when 
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r=3.68 for model P (a steep spectrum) im- 
plies a significant thermal plasma contribu- 
tion, but model M again implies an unusu- 
ally low abundance of Z = 0.003 Zq. As for 
the halo fits, we argue that such a low ap- 
parent abundance can result from a model 
that is too simple. To find a more physi- 
cally reasonable description of the data we 
examine more complex models choosing a 
priori values of Z = 1.0 Zq and Z = 0.5 
Zq for the abundance. Both are consistent 
with the data. 

For solar and one-half solar abundances, 
the data are well described with a model 
that consists of a power-law and MEKAL 
component (Model PMl, Table 1). For this 
model the best-fitting power-law index is 
r = 0.77; however when the ASCA data 
are taken into consideration, if F were re- 
ally this small it would seriously overpre- 
dict the flux above 2 keV when extrapo- 
lated to higher energies (§ 4.2). Also, if the 
"flat" hard component were thermal, the 
required temperature would be much too 
high and also inconsistent with the ASCA 
data. Alternatively, we suggest that the 
small index is an artifact of the effect of 
absorption on an intrinsically steeper hard 
X-ray component. If we allow for some ad- 
ditional absorption and fix F at the best- 
fitting ASCA value of 1.9 (Aiodel PM2, Ta- 
ble 1), we derive a column density of ~ 10^^ 
cm~^. This model is more reasonable than 
the model with the unusually small power- 
law index and it is also consistent with the 
images since we see the absorption band 
due to the galaxy disk cutting in at around 
0.7 keV. 



field background is subtracted, and even tiiougli we do 
not know tiie geometry of the system or liow much of the 
diffuse halo emission is projected onto the disk a priori, 
the fits using local background are more reasonable and 
better match the imaging data. Thus we rely on the latter 
method. 



To confirm the significant absorption in the 
galaxy disk we examine the northern disk 
(ND) and southern disk (SD) separately. 
Both spectra are plotted in Figure 4. Sub- 
tracting field background, the spectra have 
a similar shape even though there is a clear 
band of absorption in the image, which con- 
firms our suspicion that diffuse halo emis- 
sion is projected onto the band of underly- 
ing absorption in the disk. When the local 
(diffuse) background is subtracted A^h(ND) 
is 2 to 3 times A?h(SD), which is more in 
line with the image. Depending on the 
choice of models (power law or MEKAL), 
the intrinsic absorption for the hard com- 
ponent, i.e., due to the disk, can be as large 
as 6 X 10^^ cm-2. 

3. The Core Emission: 

Single-component power-law or MEKAL 

models provide poor fits to the core spec- 
trum (models P and M, Table 1). On 
the other hand, a two-component MEKAL 
plasma {Z ^ 0.5 Zq) plus power-law model 
provides an excellent fit with x^/^ — 27.6/24 
(Table 1, model PMl). In this case, F(PSPC) 
is smaller than r(ASCA), a result similar to 
the hard model component of the disk emis- 
sion. The small index again implies signif- 
icant absorption. Fixing F at the ASCA 
value of 1.9 (model PM2) yields iVn = 5 
xlO^° cm~^. Constraining the hard com- 
ponent in this way also causes the fitted 
temperature of the thermal component in 
the PSPC to be equal to the temperature 
derived from ASCA. This provides a nice 
consistency check and allows us to identify 
the soft component in the core spectrum 
with the warm ASCA thermal component 
(P97). If the hard component is more ab- 
sorbed than the thermal component (model 
PM3, our best fit), then A^H(core) can be 
as high as 4 to 5 X 10^^ cm~^. This con- 
firms the amount of absorption intrinsic to 
the galaxy disk and is consistent with that 
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inferred by VP98. 

4. Compact Sources: 

The compact sources in NGC253 are dis- 
cussed at length by VP98, who perform 
a sensitive search for point sources in the 
disk, resulting in 27 detections in the HRI 
image. We do not perform a similarly de- 
tailed analysis because our intent is only 
to identify the major point-source contri- 
butions in the PSPC data (Fig. 3 and pa- 
per 1). Our integrated PSPC point-source 
luminosity is ~ 2 x 10^^ ergs s^^, which 
is comparable to that of VP98, so wc are 
confident that we have identified the ma- 
jority of the flux from the point sources. 
Anything else is not important for our in- 
vestigation. 

The "hard" X-ray point sources can be 
modeled with a power law having P ~ 2 
(model P, Table 1) or a bremsstrahlung 
spectrum with kT ~ 1.9 (model B). The in- 
dex of the power law is similar to T{ASCA) 
for the hard component and implies that 
the i?osai-identified point sources comprise 
a significant portion of the hard X-ray por- 
tion of the ASCA spectrum. Our result 
suggests that much of the hard X-ray emis- 
sion is truly point-like, in contrast to the 
claim of C99 that the 3-10 keV emis- 
sion is primarily from hot, diffuse gas in 
the galaxy disk. From extrapolation, the 
compact sources, including the core, pro- 
vide 30 - 50% of the total 2-10 keV emis- 
sion seen with ASCA and so must make up 
a significant fraction of the 3 — 10 keV "ex- 
tent" seen with BeppoSAX (C99). 

The "soft" X-ray sources that are point-like 
in the HRI and the PSPC are located ex- 
clusively in the halo and have spectral char- 
acteristics resembling the diffuse gas (Ta- 
ble 1). This suggests that the point-like 
regions in the halo are not distinct com- 
pact sources but might be areas of locally 



enhanced halo emission, as also suggested 
by VP98. For the same two-temperature 
plasma model applied to the halo emission 
(model MM), the cool and warm compo- 
nents contribute about 40% and 60% of 
the flux, respectively, similar to the ratio 
of these components in the halo. 

We conclude from our analysis of the spatially- 
resolved PSPC spectra that there is significant 
absorption within the galaxy disk of NGC 253 on 
the order of 1 to a few x 10^^ cm~^ (in agreement 
with VP98). We find the PSPC point sources to 
be a significant source of X-rays at 2 keV and 
they are likely to contribute 30 — 50% of the total 
flux in the 2 — 10 keV band, i.e., in the ASCA 
and BeppoSAX spectra and the BeppoSAX image 
reported by C99. 

4.2. Comparing PSPC and ASCA Spec- 
tral Modeling 

Here we discuss the method of spatial-spectral 
iteration by which we derived a self-consistent 
modeling of the ASCA and PSPC spectra. In 
the above analysis, we used simple models to flt 
the spatially resolved PSPC spectra in our at- 
tempt to infer the spatial origins of the ASCA 
emission. The basic assumption is that if a PSPC 
spatial component is mostly hard or mostly soft, 
then this region can be associated with the cor- 
responding hard or soft component in the ASCA 
spectrum. Conversely, if we do not find a PSPC 
spatial/spectral component that "matches" one 
of the ASCA spectral components, then we know 
that at least one of our model assumptions is 
wrong. 

The ASCA spectrum clearly contains strong 
line emission from O, Ne, Fe, Mg, Si, S, and Ar 
(P97 and others). The integral spectrum can be 
well fitted with the sum of an absorbed hard X- 
ray component (power law or bremsstrahlung) 
and one or more thermal plasma components. 
But this is not necessarily the correct model. 
The data are ambiguous in that the amount of 
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absorption for the various components and the 
abundance of the gas can differ by orders of mag- 
nitude depending on the input assumptions. 

Figure 5 is a guide to show which component 
of the X-ray emission, hard or soft, dominates the 
spatial components in the PSPC and hence which 
ASCA and PSPC models yield consistent results. 
Since our goal is to match up ASCA and the 
PSPC, we choose the two most significant ASCA 
components from P97 and compare their spectral 
parameters with the PSPC. Figure 5a shows F vs. 
A^H for the power-law (hard X-ray) ASCA com- 
ponent compared with the power-law results for 
the spatially-resolved PSPC spectra. All PSPC 
points are derived for a single power-law model 
except for core(2), which represents the photon 
index for the best-fitting two-component core 
model (model PM3, Tabic 1). The hard X-ray 
point sources and the hard component in the 
PSPC core clearly have the spectra that are most 
similar to ASCA. This consistency provides inde- 
pendent support for our complex modeling of the 
core spectrum. 

Figure 5b shows kT vs. Nu for a MEKAL fit 
to the ASCA soft component (P97) and single- 
component MEKAL fits to the PSPC, except 
for core(2) as above. The disk emission with 
kT ~ 0.8 keV most closely corresponds to the 
soft ASCA thermal component found by P97; 
however, in this case, no spatially resolved PSPC 
component has exactly the same model parame- 
ters. In fact, much of the spatially resolved emis- 
sion has a significantly lower temperature than 
measured with ASCA. These results strongly im- 
ply a three-component modeling of the PSPC 
and ASCA data, including a very soft ther- 
mal component that contributes in the ASCA 
band. For comparison we also plot the single- 
component model result for the integral PSPC 
spectrum, which illustrates the incorrect results 
we might infer if the data had poor enough statis- 
tics that a single-component model were statisti- 
cally acceptable. 



4.3. Modeling the ASCA Spectrum 

We compare our results for two- and three- 
component models for ASCA in Table 2. To ac- 
count for systematic differences between our fit- 
ting techniques and those of P97, we have re- 
extracted the ASCA spectrum using their re- 
gion sizes. The abundances of individual ele- 
ments are not well determined so we follow P97 
and group the plasma abundances as follows: 
Ne and Fe (hereafter Z-^-p), O, Mg, Si, and S 
(hereafter ^oMSs); and N, Ar, Ca, and Ni (here- 
after ^nacn)- In reality, Ne should be tied to 
S, Si, and Mg, based on the physics of produc- 
tion in Type II supernovae; however, the fact 
that the Ne emission is blended with Fe L (Fig- 
ure 1), combined with the poor energy resolution 
of the detectors, forces us to tie this parameter 
to Fe. This leaves 13 free parameters for the 
ASCA fits, which arc the normalization of the 
cool component^, the temperature (kT(w)), nor- 
malization, .^oMSS) -^NF) -^NACN) and the absorb- 
ing column (iVH(w)) of the warm component, the 
index (or temperature), normalization, and ab- 
sorbing column (A^H(h)) of the hard (or "hot") 
component, and the relative normalizations of 
the instruments. The Galactic foreground col- 
umn density of 0.9x10^^ cm~^ is included in 
all cases. Based on joint PSPC and ASCA fits 
(§ 4.4) we choose Znf,-^nacn = 1-0 Zq and fix 
the temperature and abundance of the cool gas 
at kr(c) = 0.2 keV and Z = 0.2 Zq, respectively. 

We derive results similar to P97 for a two- 
component model (models PM or BM in Table 2) 
but find equally good or better fits for a three- 
component model (models PMM or BMM). We 
conclude that the three-component and two-component 
models are statistically equivalent (the intrinsic 
uncertainty in xLin results from the choice of 
plasma models) and ambiguous. But when we 

''We do not measure multiple abundances for the cool com- 
ponent because the line emission for the ~0.2 keV plasma 

is dominated by Ne and Fc L and the lines at higher en- 
ergies are swamped by the emission from the warmer gas. 
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account for the low-temperature component in 
the PSPC, the three-component model is much 
more consistent with the imaging data. The 
three-component model not only allows the metal 
abundances of the warm component to be 1 — 2 
^0 (which is physically more meaningful for a 
galaxy with a very high star formation rate, cf. 
Zaritsky, Kennicutt and Huchra 1994), it also al- 
lows A^H(h) to be smaller and more consistent 
with A^H inferred from PSPC fits. Our best- 
fitting ASCA model (PMM, Table 2) is shown 
in Figure 6. 

4.4. Joint Spectral Modeling 

Having demonstrated the large uncertainty in 
deriving absolute abundances for thermal mod- 
els, we fix Znacn at values of 0.5 or 1.0 ^Qfor 
the joint ASCA and PSPC fits. With ^nacn = 
0.5 Zq (Table 3), the best fit is obtained with 
two-temperature MEKAL plasma plus a power 
law (x^ = 538.5 for 499 d.o.f., model PMMl). 
For this model we derive .Z^OMSS = 0-8 and Znf 
= 0.4. However, space is shallow in the region 
near Xmin due to Znf, -^omss, Nii{h), and iVH(m) 
being strongly coupled. After thorough testing of 
the allowed parameter space, we conclude that 
there is no need for the Fe abundances to be 
less than solar. The joint data and best-fitting 
two-temperature MEKAL plasma plus power law 
model with Znf = 0.5 Zq are shown in Figure 7 
(model PMM2). 

In general, F, kT(c), kT(w), and kT(h) are 
well determined for the joint fits, while there are 
large uncertainties in Nn{w), Nji{h) and .^OMSS- 
ZoMSS tends to be larger than .^nFj but there 
are cases where the two are consistent to within 
the statistical errors. 

The mean value for A^h(w) is 5.81 x lO'^^ cm^^ 
and the mean value for ArH(h) is 6.72 x 10^^ 
cm~^. Both are consistent with A^H(core) in- 
ferred from the PSPC, but this is slightly larger 
than we might expect for A^h(w). If we add 
a third MEKAL plasma component with kT = 
0.14 keV to match the softest emission in the 



PSPC, then N}i{w) can be lower and more con- 
sistent with the PSPC data. This also allows 
the abundance in the halo to be higher than 0.2 
Zq. An alternative explanation of an effectively 
large Nn (w) is a more complex absorption struc- 
ture. If the hard X-ray component is modeled as 
being partially covered, about 30% of the hard 
X-ray flux could be unabsorbed or absorbed by 
a very low column [model (PMM)pc. Table 3]. 
In other words, if some of the hard X-rays are 
passing through a patchy absorber, then it would 
add flux at low energies, another reason why we 
might measure low abundances below 1 kcV for 
purely thermal models for the soft emission. 

An independent estimate of NhQli) and N}i{w) 
arises from comparing the PSPC and ASCA 
fluxes. The mean fluxes for the joint model, in 
units erg s~^ cm~^, are f(c) = 2.88 xlO~^^, f(w) 
= 1.03 X 10-12^ and f(h) = 1.13 xlO^^^ ^an 
compare these with the PSPC fluxes of f(c) = 
~ 2.0 X 10-12^ f(w) = 1.72 xlO-i^ and f(h) = 
1.91 xlO~^^ erg s~^ cm~^. The column density 
that produces the best agreement between the 
fluxes of the hard and warm components is 4 to 
5 xlO^^ cm~^, which is reasonable for a highly- 
inclined spiral galaxy disk. 

4.5. Summary 

Our conclusions for NGC 253 are: 

(1) There is a large amount of intrinsic absorp- 
tion in the disk of the galaxy {Nu^ 1 — 6 x 10^^ 
cm~^). This absorption complicates spectral flts. 

(2) The abundance measurements from ASCA 
and PSPC spectra are highly uncertain but are 
statistically consistent with solar values. Also, 
the a-burning elements are not required to be 
unusually overabundant relative to Fc. 

(3) A three-component model that consists 
of a power law and two MEKAL plasma com- 
ponents with solar or near solar abundances 
provides a fit as good or better than a two- 
component model consisting of a power law and 
a MEKAL plasma having unusual abundances. 
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To maintain consistency with spatially-resolved 
emission, the first alternative is allowed by the 
Rosat data; the second is not. 

(4) Extrapolating the PSPC results into the 
ASCA band allows us to infer which spatial com- 
ponents dominate the ASCA spectrum. The 
compact sources, including the core, provide 30 — 
50% of the total 2 — 10 keV emission seen with 
ASCA and so must make up a significant fraction 
of the 3 - 10 keV "extent" seen with BeppoSAX 
(C99) when the telescope point-spread function 
is accounted for. 

Comparing our results with those by C99 it 
is important to realize that our three-component 
model is a composition based on emission mod- 
els describing X-ray emitting objects known to 
exist in virtually all spiral galaxies. While C99 
invoke the presence of an additional, extremely 
hot thermal component, they neglect all poten- 
tial continuum emission in the 2 — 10 keV band 
from XRBs. However, there is no reason to be- 
lieve that NGC253, with its massive starburst, 
should not host any X-ray emitting massive bi- 
naries. On the other hand, XRBs can account 
for the observed thermal Fe K line emission. 

5. Spectroscopic Results for M82 
5.1. Modeling the PSPC Spectra 

1. The Diffuse Halo Gas: 

Power law or single-temperature thermal 
plasma models are statistically ruled out 
for the halo emission in M 82 (models P and 
M, Table 5). Assuming solar abundances, 
the halo emission is best described with 
a two-temperature model having kTi = 
0.3lt[!:!]5 keV and kTa = 4..70lfs2 kcV, re- 
sulting in = 14.3 for 25 d.o.f. (model 
MM, Table 5). As in the case of NGC 253, 
we cannot rule out a power law as the 
harder of the two components instead of a 
thermal plasma, although the thermal de- 
scription is somewhat preferred from the 
ASCA result (§ 5.2). 



2. The Core Spectrum: 

Two model components are required to fit 

the core spectrum. If these are a power law 
and a MEKAL plasma (model PMl), then 
the hard component prefers to be much 
flatter than measured with ASCA (simi- 
lar to NGC 253). This could result from 
significant unresolved line emission above 
1 keV (Si) or from a heavily absorbed 
component. When the hard component is 
constrained to have the same photon in- 
dex as the ASCA hard component (model 
PM2, identical to the technique we used for 
NGC 253), the inferred thermal tempera- 
ture in the core approaches that in the halo. 
In this case, the hard component requires 
absorption in excess of the Galactic value. 

3. The Compact Sources: 

Most of the point-like sources near M 82 do 
not appear to be associated with the gal- 
axy since they fall outside the plane of the 
galaxy. However, the very soft sources to 
the north have properties similar to the hot 
halo gas in M82 and they may, in fact, 
be associated with the outflow (Lehnert, 
Heckman and Weaver 1999). See paper 1 
for further discussion of the point sources. 

5.2. Modeling the ASCA Spectrum 

ASCA results for M 82 have been published by 
P97, T97, and ML96. All are of comparable sta- 
tistical quality, but we propose that the spatial 
information from the PSPC can be used to tell 
which are physically meaningful. 

Similar to NGC 253, we performed a detailed 
investigation of the ASCA spectrum using a 
two-component "baseline" model to determine 
X^(min) and comparing this with the three-component 
model. Since the statistics are good enough, 
we also use both Raymond-Smith and MEKAL 
plasma models to see how well the data can dis- 
tinguish between the two. For a model that con- 
sists of a thermal plasma with solar abundances 
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and a power law, all fits are poor with Xv 
ues ranging from 1.44 to 1.48 (Table 6, models 
PRs and PMg). We next allow the abundances 
of individual elements to vary, although keeping 
N and Na at their solar values because of their 
large uncertainty. In this case, the fits improve 
dramatically (models PRy, PMy and BMy) and 
we derive abundances similar to P97. 

Like NGC 253, the PSPC data clearly indicate 
an additional low-temperature plasma and excess 
absorption (Strickland ct al. 1997). However, 
three-component models with solar abundances 
(models PRRg, PMMg) do not fit the ASCA data 
as well as two-component models with variable 
abundances. In fact, when the abundances are 
allowed to vary, a soft component is not required 
in the ASCA bandpass at all. 

We also detect Fe K emission in M82, which 
could have a thermal origin or could come from 
X-ray binaries. Adding a narrow Gaussian to 
the ASCA model yields Ax^ = 11 and a line 
at 6.65 ± 0.10 keV with an equivalent width of 
77±43 eV and normalization of 1.63±0.90x 10"^ 
photons cm~^ s~^. 

5.3. Joint Spectral Modeling 

For joint ASCA and PSPC fits of M 82, we are 
forced to ignore the ASCA data below 0.7 keV to 
avoid a slight miscalibration between the detec- 
tors. Currently, there is a calibration effect in 
the SIS at energies less than 1 keV that causes 
an apparent decrease in the quantum efficiency 
and mimics excess absorption on the order of a 
few xlO^° cm~^. The magnitude of the problem 
depends on the SIS CCD mode and the time of 
the observation since launch (Weaver and Gel- 
bord 1999; in prep.). Early in the mission, this 
effect is only detectable for fairly bright sources 
and so is not a problem for our joint analysis of 
ASCA and PSPC spectra for NGC 253 (§ 4.4). 

In the case of M 82 our best fits require the ad- 
dition of a second low-temperature plasma com- 
ponent (models PMMg and PMMy, Table 8), re- 



gardless of whether or not the abundances are 
allowed to vary. For variable abundances, O and 
Mg are uncertain by a factor of ~5 and the ratio 
of O to Fe ranges from 0.92 to 3.6. Non-solar 
abundances are required for some elements, but 
we find that the abundances are not necessar- 
ily unusual or overly large when excess absorp- 
tion is allowed. This result differs from T97 and 
P97, who instead derive or assume a low value 
for N}i{w). Our analysis clearly suggests that 
the relative abundances are model dependent and 
difficult to determine. 

The most appropriate self-consistent descrip- 
tion for the PSPC and ASCA data is a three com- 
ponent model that consists of cool thermal emis- 
sion, warm thermal emission and a hard compo- 
nent. Model PMMy is shown in Figure 8. 

6. Summary &; Implications 

6.1. Comparison to Previous Work 

We have arrived at somewhat different con- 
clusions than have several other recent analyses 
of the X-ray spectra of starbursts. As we have 
argued above, this is not surprising. Starburst 
galaxies are complex X-ray sources, with emis- 
sion resulting from diffuse, hot gas inside and 
outside the galaxy (associated with starburst- 
driven winds). X-ray binaries (XRBs) and super- 
novae in the disk of the galaxy, and absorption 
from the galaxy. The soft X-ray obscuration, in 
particular when combined with limited spatial in- 
formation, causes ambiguity in unraveling the de- 
tailed physics of the hot gas (e.g., abundances), 
because a multi-temperature distribution of soft 
X-ray thermal spectra can be modified by under- 
lying absorption and can imitate the signature of 
gas with unusual abundance ratios. 

From ASCA and Rosat PSPC archival data for 
NGC 253 and M82, we find ambiguities in the 
spectral model fitting such that multiple com- 
ponents are confused when there is no adequate 
spatial information. Our most plausible model 
that is consistent with the PSPC images and 
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spectra consists of three components, (1) a low- 
temperature thermal plasma, (2) a warm ther- 
mal plasma, and (3) a hard power-law or thermal 
component. We identify the spectral components 
with the likely emitting sources in Table 9. The 
line-of-sight absorption is approximately a few 
times 10^^ cm~^ in the disk of the galaxy, but 
both this and the plasma abundances are model- 
dependent, with most constraints being inferred 
from the imaging data. Our self-consistent ex- 
amination illustrates the criticality of the X-ray 
images when searching for a physically reason- 
able description of the data. The absolute abun- 
dances are impossible to determine without the 
improved spatial/spectral resolution of Chandra 
and XMM. 

Our results differ from those by other authors. 
T97 examined a three component model, but 
fixed the absorption of the thermal components 
to be significantly lower than measured in the 
disk (VP98). They also choose to investigate 
abundances and non-equilibrium effects rather 
than absorption effects, which we claim have a 
dramatic impact on the soft X-ray spectrum. 
While none of the spectral models used by us, 
P97, T97, or others are unique interpretations of 
the data, it is worth emphasizing that there is 
no overwhelming evidence at this time to support 
highly unusual abundances or non- equilibrium ef- 
fects for M82 or NGC253. 

Joint ASCA + PSPC fits prefer larger absorb- 
ing column densities than the PSPC alone, ex- 
cept perhaps with a bremsstrahlung modeling of 
the hard component. This appears to be an ef- 
fect of the different ways in which the fits are 
weighted, with the PSPC fits weighted toward 
the soft emission that is not heavily absorbed 
and the joint fits weighted more toward the warm 
emission that is absorbed. We expect the joint 
fits to be more sensitive to intrinsic, underly- 
ing absorption of the order 10^^ cm~^ than the 
PSPC, because such an absorbed component is 
hard to isolate in the PSPC data. 

The ASCA data are barely of sufficient qual- 



ity to comment on the spectral origin of the hard 
X-ray component, i.e., thermal or non-thermal. 
We detect weak Fe K emission in M82, which 
could arise from thermal processes or from XRBs. 
There are disagreements about a thermal inter- 
pretation for the hard X-rays in NGC 253 be- 
cause although BeppoSAX has now detected the 
expected Fe K line in NGC 253, the emission 
is weak (P97 and Ohashi ct al. 1990). One 
explanation for the weak Fe lines is that the 
hard emission contains a significant contribution 
from point sources - XRBs or SNe. Our spec- 
tral/spatial analysis of NGC 253 suggests that 
compact sources must at least partly contribute 
to the hard emission. 

6.2. The Origin of the Thermal X-ray 
Emission in Starbursts 

It has long been known that starbursts are un- 
usually intense X-ray sources (cf., Fabbiano 1989 
and references therein), and likewise, the exis- 
tence of X-ray-emitting superwinds is now well- 
established. However, the process (es) by which 
they produce X-ray emission is (are) not clear 
(see, for example Strickland 1998). The chem- 
ical abundances in the X-ray-emitting gas can 
provide strong constraints on these processes. 

If the outflowing wind consists purely of the 
thermalized ejecta from supernovae and stellar 
winds, the wind fluid would have extremely high 
metallicity (~ 10 Zq), with the ratio of a-burning 
elements to Fc roughly twice solar (e.g. Gibson, 
Loewenstein, & Mushotzky 1997). In this case, 
the gas would be very hot (~ 10^ K prior to adi- 
abatic cooling), but so tenuous that it will be an 
insignificant source of X-ray emission (Chevalier 
& Clegg 1985; Suchkov et al 1994). Thus, the 
observational constraints on the absolute metal 
abundances in the X-ray plasma agree with other 
arguments: we arc not observing emission from 
undiluted type II supernova ejecta. 

Various processes have been proposed to boost 
the X-ray luminosity produced by a superwind. 
Substantial quantities of material could be mixed 



12 



into the wind in or near tlic starburst as the stel- 
lar ejecta interact with the starburst ISM (the 
wind could be centrally "mass-loaded" - Suchkov 
et al 1996; Hartquist et al 1997). In this case 
the abundances will be intermediate between the 
extreme values for the supernova ejecta and the 
more normal values for the gas in/near the star- 
burst. For example, the hot thermal compo- 
nent we observe in the cores of NGC253 and 
M82 has a temperature of roughly lO'^ K, which 
could imply a roughly 10:1 dilution of the hot 
(10* K) supernova ejecta with ambient gas. In 
this case, if the ambient gas had solar abun- 
dances, the X-ray gas would have abundances of 
a-elements like O, Ne, Si, Mg that are roughly 
twice solar and only slightly super-solar a/Fe ra- 
tios. Such abundances are consistent with our 
three-component fits to the joint X-ray spectra 
of M 82 and NGC 253. 

Alternatively, very little central mass-loading 
might occur, and the hot, tenuous (X-ray faint) 
wind fluid could propagate into the ISM and 
galactic halo where it could shock heat the halo 
material to 10*^ to lO'^ K, causing it to emit 
soft X-rays (e.g., Suchkov et al 1994). In this 
case we would expect the abundances in the X- 
ray-emitting gas to reflect the pre-existing abun- 
dances in the ambient gas. In the disk of the 
"host" galaxy or within the starburst itself, this 
model would predict roughly solar abundances 
(consistent with our three component spectral 
models). In the halo, substantially subsolar 
metal abundances are possible. The High Ve- 
locity Clouds in the halo of the Milky Way gen- 
erally have metal abundances less than 10% so- 
lar (Wakker &: van Woerden 1997). The Lyman- 
Limit quasar absorption-line systems are gener- 
ally interpreted as arising in the halos of nor- 
mal galaxies (Steidel 1993), and these have abun- 
dances ranging from < 0.01 to 0.4 solar (Prochaska 
1999; Prochaska & Buries 1999). Thus, i/ the 
very low metallicities implied by the simple two- 
component fits are correct, it would imply that 
the observed X-ray emission is arising from shock- 



heated halo gas. This may be tenable for the dif- 
fuse halo emission in NGC 253 and M 82, but not 
for the emission in the starburst core. 

It is also worth emphasizing that we are now 
starting to uncover a complexity in the ISM of 
starbursts that is similar to that in the Milky 
Way. This increasing complexity (e.g., two or 
three thermal components in the hot X-ray-emitting 
gas alone) supports the suggestions by Norman & 
Ferrara (1996) that a more general description of 
the ISM (a phase-continuum) should be used in- 
stead of simply adding more and more individual 
phases. Strickland (1998) has advocated this ap- 
proach in fitting the X-ray spectra of starbursts. 

6.3. Starbursts & the Enrichment of the 
Intergalactic Medium 

The results in this paper imply that - con- 
trary to earlier analyses - the hot X-ray-emitting 
superwind gas can have normal (roughly solar) 
abundances, and therefore may be carrying a sig- 
nificant amount of metals out of the starburst 
galaxy and into the intergalactic medium (IGM) . 
This corroborates the suggestion by Heckman, 
Armus, k Miley (1990; hereafter HAM90) that 
starburst-driven (i.e., mostly dominated by type 
II supernovae - hereafter SN II) outflows are one 
of the potential processes for enriching the IGM 
at the present epoch. The importance of this ef- 
fect in comparison with others is currently a mat- 
ter of debate. While the observed values of the 
a/Fe ratios in the intra-cluster medium (ICM) 
in galaxy clusters leave no doubt that most of 
the metals (by mass) were contributed by SN II 
(Gibson, Loewenstein, & Mushotzky 1997; Ren- 
zini 1997), the roles of Type I supernovae, galac- 
tic winds, tidal and ram-pressure stripping, ellip- 
tical vs. spirals, and dwarfs vs. massive galaxies 
are still matters for debate. 

For example, Gnedin (1998) argues that tidal- 
stripping of the ISM is the dominant mechanism 
by which the IGM was "polluted", and that the 
role of SN Il-dominated outflows is negligible on 
cosmological timescales. On the other hand, oth- 
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ers have argued that bimodel star-formation in 
elhptical galaxies, i.e., high-mass star- formation 
with an SN Il-dominated ISM in the early phase 
that induces a massive wind and dies after about 
0.04 Gyr (Arnaud et al. 1992; Elbaz et al 1995) 
might have caused significant enrichment of the 
IGM from ellipticals then, with no observational 
trace of further significant star-formation at the 
present epoch (see Zepf & Silk 1996). 

Wiebe et al (1999), on the other hand, state 
that disk galaxies contribute significantly to the 
enrichment of the IGM, with an "'effective" loss 
(per unit luminosity) comparable to that of el- 
lipticals, which dominate the enrichment of the 
IGM only because they outnumber spirals. The 
efficiency of metal losses of spirals into the IGM 
depends on their total mass. Nath k. Chiba 
(1995) find that dwarf starbursts can contribute 
to the metal enrichment of clusters at high z, but 
they produce only low-metallicity gas. 

Kauffman &; Chariot (1998) have proposed 
that hot winds of merging galaxies at high z, 
which later form giant ellipticals, are the sources 
of metals in clusters. This idea could logically 
combine the two apparently contradicting models 
favoring metal ingestion into the IGM by either 
ellipticals or spirals. The question remains how 
significant the metal enrichment of the IGM due 
to galactic superwinds in spirals is at the present 
epoch. 

The most likely analogs to these systems in the 
current universe are far infra-red luminous galax- 
ies (FIRGs), like e.g. NGC6240 and Arp220 
(HAM90). According to these authors, FIRGs 
could inject of order 10^ Mq and 10^^ ergs per 
-L* galaxy over a Hubble time, with no evolu- 
tion in the superwind rate over time. While 
these systems are very scarce at low redshifts, 
sub-millimeter observations show that dusty ul- 
traluminous galaxies are far more common at 
high-redshift, and may rival normal star- forming 
galaxies in terms of metal production over the 
age of the Tinivcrsc (e.g., Lilly 1999). More- 
over, galaxies like M82 and NGC253 exist in 



far greater numbers than ultraluminous mergers 
in the current universe, and could dominate the 
present-day enrichment of the IGM in groups and 
clusters. 

Material torn out of interacting galaxies can 
be observed in the form of tidal tails and sim- 
ilar structures (Toomre & Toomre 1972). The 
existence of such structures is a tracer of the 
shallowness of the galaxies' gravitational poten- 
tials, which also makes it easier for them to ex- 
pel metal-enriched gas via supernova-driven su- 
perwinds (see, e.g., NGC4631; Weliachew et al 
1978, Combes 1978, paper 1). In this context it is 
important to understand the abundances in these 
outflows. We have found in paper 1 that super- 
winds, with hot gas up to 10 kpc above the disks 
of galaxies, do exist and that outflows are a com- 
mon phenomenon in nearby starburst galaxies. 
The current paper shows that the metallicities of 
the outflowing gas can be high enough to allow 
for significant metal losses into the IGM, if escape 
velocity is reached (cf., estimates by HAM90). 

The observational limit is currently that the 
existing ambiguity in the minimum space of 
the X-ray spectral fits cannot be removed. At 
present we can only show that it is indeed pos- 
sible to fit the data of both M 82 and NGC 253, 
two of the IR-brightest galaxies in the nearby 
Universe, in accordance with chemical evolution- 
ary models and that this implies possible signifi- 
cant metal losses into the IGM. However, at the 
same time spectral fits to the data with extremely 
low metallicities are still allowed. As the next 
step, data of better quality is needed, first for 
these nearby galaxies and then for a larger sam- 
ple of galaxies to see whether one of the two solu- 
tions can then be excluded. This will be possible 
with the next generation of X-ray telescopes, like 
Chandra or XMM. 

X-ray studies of distant (high-z) starbursts 
can be used to estimate the IGM enrichment 
by starburst-driven superwinds as a function of 
epoch. However, this will also only be possi- 
ble based on data from one of the future X-ray 
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missions. A good diagnostic tool is the [0/Fc] 
abundance ratio, which should be higher than 
Solar for SN Il-dominated winds, while for an 
SN I-dominated ISM it should be close to Solar 
(e.g., Matteucci & Vettolani 1988). Future high 
throughput X-ray observations should put much 
tighter constraints on the possible metallicities of 
the hot gas in star-forming galaxies and thereby 
help to shed light on the potential metal enrich- 
ment of the IGM as a consequence of massive 
star-formation. 

6.4. Current Limitations &; Future Prospects 

Combining the data from Rosat PSPC and 
ASCA is a step forward in our attempt to de- 
scribe the complex X-ray emission properties of 
galaxies like NGC253 and M82 compared to 
studies based on data from only one instrument. 
Not only is there an ambiguity in the choice of the 
"best" model for the combined spectra, but an 
interpretation of data from only one satellite (in 
the pre- Chandra and pre- XMM era) can be even 
more misleading because different bandpasscs are 
sensitive to different emission components, and 
neither Rosat nor ASCA can detect all of them. 

X-ray spectral fits now turn into multi-component 
multi-d.o.f. approximations including internal (X- 
ray imaging) and external (observations in other 
wavebands) consistency checks in the choice of 
model components, as performed here. At present, 
the choice of models to be fitted to the data 
depends on assumptions or information from 
sources other than the spectra themselves. 

However, there are still a few shortcomings, 
which call for the improved capabilities of the 
next generation of X-ray satellites: 

• While Rosat has good spatial resolution 
but limited spectral resolution, ASCA has 
a poor spatial resolution but high spec- 
tral resolution. This leads to possible am- 
biguities and inconsistencies in combining 
the data. The superb spatial resolution of 
Chandra, coupled with an energy resolu- 



tion similar to ASCA should allow the spec- 
tral properties to be far better determined, 
since it will be possible to spatially segre- 
gate the different emitting components. 

• The effective areas of the past generations 
of X-ray telescopes were not high enough to 
provide adequate signal-to-noise for the in- 
dividual X-ray sources/components in star- 
bursts (especially given that M 82 and NGC 253 
are at least an order-of-magnitude brighter 
X-ray sources than any other starburst gal- 
axy). XMMs combination of good angular 
resolution, large effective area, and good 
energy resolution will be a major improve- 
ment. 

• For gas in the typical temperature range 

seen in starbursts, even the energy resolu- 
tion of the CCD's on XMM and Chandra 
are inadequate to resolve the Fe-L complex 
(whose strength and mean energy largely 
drive the fits for temperature and metal 
abundance in the existing spectra). The 
~ 10 eV energy resolution of the XRS in- 
strument on Astro- E should allow many in- 
dividual X-ray lines to be resolved. For 
sources as bright as M82 and NCC 253, it 
may be possible to measure the tempera- 
ture range and chemical abundances via a 
multitude of independent emission-line ra- 
tios (as is done in solar /stellar coronae) 
rather than by fitting idealized models to 
lower-resolution, complex spectra. 

• Both Chandra and XMM cover a very wide 
energy range that includes the bandpasses 
of Rosat + ASCA, thus removing the un- 
certainties connected with the combination 
of two datasets with disparate properties. 

Acknowledgements We thank the referee, Joel 
Bregman, for helpful comments. This research 
was supported in part by NASA grants NAG5- 
6400 and NAG5-6917. 



15 



REFERENCES 

Arnaud, M., Rothenflug, R., Boulade, O., Vi- 
groiix, L., & Vanioni-Flam, E. 1992, A&A, 
254, 49 

Bregman, J. N., Schulman, E., &; Tomisaka, K. 
1995, ApJ, 439, 155 

Cappi, M. et al. 1999, A&A, 350, 777 

Chevalier, R. & Clegg, A. 1985, Nature, 317, 44 

Combes, F. 1978, A&A, 65, 47 

Dahlem, M., Weaver, K. & Heckman, T. 1998, 
ApJS, 118, 401 (paper 1) 

Elbaz, D., Arnaud, M., Casse, M. et al 1995, 
A&A, 265, L29 

Fabbiano, G. 1982, ApJ, 330, 672 

Fabbiano, G. 1989, ARAA, 27, 37 

Fabbiano, G. & Trinchieri G. 1984, ApJ, 286, 491 

Gibson, B., Loewenstein, M., &; Mustiotzky, R. 
1997, MNRAS, 290, 623 

Gnedin, N. Y. 1998, MNRAS, 294, 407 

Hartquist, T., Dyson, J., & Williams, R. 1997, 
ApJ, 482, 182 

Heckman, T. M., Armus, L., & Miley, G. K. 1990, 
ApJS, 74, 833 (HAM90) 

Kaastra, J.S. 1992, An X-Ray Spectral Code 
for Optically Thin Plasmas (Internal SRON- 
Leiden Report, updated version 2.0) 

KaufFman, G., & Chariot, S. 1998, MNRAS, 294, 
705 

Lehnert, M. D., Heckman, T. M., & Weaver, K. 
A. 1999, ApJ(1999astro.ph..4227L) 

Lilly, S. J. et al. 1999, ApJ, 518, 641 

Matteucci, F., & Vettolani, G. 1988, A&A, 202, 
21 



Mcwc, R., Gronenschild, E.H.B.M., and van den 
Oord, G.H.J. 1985, A&AS, 62, 197 

Mewe, R., Lemcn, J.R., and van den Oord, 
G.H.J. 1986, A&AS, 65, 511 

Moran, E. C, & Lehnert, M. D. 1997, ApJ, 478, 
172 

Nath, B. B., & Chiba, M. 1995, ApJ, 454, 604 

Norman, C. A., & Ferrara, A. 1996, ApJ, 467, 
280 

Ohashi, T., Makishima, K., Tsuru, T., Takano, 
S., Koyama, K. & Stewart, G. C. 1990, ApJ, 
365, 180 

Petre, R. 1993, BBXRT: A Preview to Astronom- 
ical X-ray Spectroscopy in the 90' s (NASA: 
Goddard Space Flight Center), p. 321 

Prochaska, J. 1999, ApJ, 511, L71 

Prochaska, J., and Buries, S. 1999, AJ, 117, 1957 

Ptak, A., Serlemitsos, P., Yaqoob, T., & 
Mushotzky, R. 1997, AJ, 113, 1286 

Raymond, J. C. & Smith, B. W. 1977, ApJS, 35, 
419 

Renzini, A. 1997, ApJ, 488, 35 

Schaaf, R., Pietsch, W., Biermann, P. L., Kron- 
berg, P. P., & Schmutzler, T. 1989, ApJ, 336, 
722 

Steidel, C. 1993, The Environment and Evolution 
of Galaxies^ ed. J.M. ShuU and H. A. Thron- 
son, Jr. (Kluwer: Dordrecht), p. 263 

Strickland, D. K. 1998, Ph.D. Dissertation, The 
University of Birmingham 

Strickland, D. K., Ponman, T. J., & Stevens, L 
R. 1997, A&A, 320, 378 

Suchkov, A. A., Balsara, D. S., Heckman, T. M., 
& Leitherer, C. 1994, ApJ, 430, 511 



16 



Suchkov, A., Berman, V., Heckman, T., k. Bal- 
sara, D. 1996, ApJ, 463, 528 

Toomre, A., & Toomre, J. 1972, ApJ, 178, 623 

Tsuru, T., Hayashi, I., Awaki, H., Koyama, K. 
& Ptak, A. 1997, RMxAC, 6, 218 

Vogler, A., & Pietsch, W. 1999, A&A, 342, 101 

Wakker, B., & van Woerden, H. 1997, ARAA, 
35, 217 

Watson, M. G., Stanger, V. & Griffiths, R. E. 
1984, ApJ, 286, 144 

Weliachew, L., Sancisi, R., & Guelin, M. 1978, 
A&A, 65, 37 

Wiebe, D. S., Shustov, B. M., k Tutukov, A. V. 
1999, A&A, (in press; astro-ph/9902202) 

Zaritsky, Kcnnicutt & Huchra 1994, ApJ, 420, 87 

Zepf, S., k Silk, J. 1996, ApJ, 466, 114 



This 2-columii preprint was prepared with the AAS 
macros v4.0. 



17 



Table 1: PSPC Spatially-Resolved Spectral Results for NGC253 



Region" 


Model^ 








Z3 








/a!(0.1-2.0keV)* 


/x(0.1-2.0kcV)'= 


















[Obs.] 


[Unabs.] 


(1) 


(2) 


(3) 


(4) 


(6) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 



Halo 
Halo 


Ml 
M2 


o-go+o-^o 

0.90+0-40 


20+0-05 

36+0-08 
0-36_o.09 


7 0+0-5 
'■"-0.6 
49 9+6.3 


1(f) 

01+0-02 
- -0.01 






96 
20, 


,1/28 
,1/27 


, ,,„+0.09 
■"-0.10 
1 3^+0.18 
1-^^-0.10 


1 85+0-13 
-°°-0.14 
„ 07+0.28 
^-"7-0.14 


Halo 


MM 


090+0.80 


14+0-03 
-0.06 

61+0-39 
- -0.39 


20.0 
5.69 


0.2(f) 
0.2(f) 






19, 


.6/26 


0.75(s) 
0.49(h) 


1.29 
0.60 


Soft src 
Soft src 


P 

M 


0.90+1-09 


40+0-12 
"-0.17 


4 2+0.7 
— 0.6 


0+0-08 
-0.0 


2 qn+0-72 
2-90_o 4g 


0.23+0.05 


14, 
11 


,1/19 

,2/lS 


0.12+0.03 

0.12+0.02 


0.31+0.06 
0.18+0.02 


Soft src 


MM 


0.99+2-50 


14+0-04 
-0.06 
0.61(f) 


1.83 
0.47 


0.2(f) 
0.2(f) 






12, 


.1/18 


0.07(s) 
0.04(h) 


0.12 
0.05 


All Disk 


P 


<i ()()+2-5 
■ - -1.8 


^ A ^+0.14 
-0.11 




r> r>r,o + 0.025 

0-003+^.^03 


, fio+0.42 
3-68_o.36 


3.7+0.3 


34, 


B/29 


1.10 + 0.10 


17.0+1.3 


All Disk 


M 


4.85+2-25 


54.5+4.1 






37, 


.5/28 


0.97+0.08 


2.46+0.19 


All Disk 


PMl 


■2.31+111 


0.26+0.06 


6.09 


0.5(f) 


0.77(bf) 


1.10 


19, 


.1/27 


0.61(s), 0.29(h) 


0.92(s), 0.33(h) 


All Disk 


PM2 


7 09+10.54 
'-"'^-3.96 


0.29+0.03 


3.83 


0.5(f) 


1.9(f) 


1.77 


29, 


,1/28 


0.48(s), 0.37(h) 


0.59(s), 0.79(h) 


N. Disk 


P 


16-514.4 








3.68(f) 


0.72+0.20 


31 


,9/25 


0.13+0.03 


3.3+0.9 


N. Disk 


M 


12.011-? 


0.44(f) 


10.4+3-3 


0.003(f) 






31, 


,1/25 


0.12+0.0 


0.47+0.13 


S. Disk 


P 


J, 24+0.63 








3.68(f) 


3.05+0.22 


24, 


.1/25 


0.96+0.06 


13.8+1.0 


S. Disk 


M 


4 33+0-52 
4-''''-0.43 


0.44(f) 


45.4+3.2 


0.003(f) 






20, 


.6/25 


0.85+0.06 


2.06+0.14 


Core 
Core 


P 

M 


^ lt;+0.67 
4 29+0-63 


, 00 + 1.17 
l-**3_o.48 


37 4+1.8 
-2.0 


09+0.17 
-0.09 


1 on+0.14 
1-80-0.12 


8.85+0.32 


108.0/27 

92.5/26 


2.02+0.06 

1 95+0-09 

-'^^-o.io 


3.67+0.13 

3 01+0-15 
''-01-0.17 


Core 


PMl 


2 T5+1-25 


0.60+0-14 
— 0.15 


2.90 


0.71{bf) 


91+0-41 
0-^-0.46 


4.99 


27, 


,6/24 


0.56(s), 1.28(h) 


0.65(s), 1.46(h) 


Core 


PM2 


"310:^1 


„ c,n+0-23 


1.64 


1.0(f) 


1.9(f) 


6.93 


70, 


.9/26 


0.33(s), 1.58(h) 


0.42(s),3.06(h) 


Core 


PM3 


2.4811:32 


51+0.12 
0-51-0.13 


13.5 


1R+0.12 
0-18-0.08 


1.9(f) 


10.6 


25, 


.2/24 


0.90(s),0. 89(h) 


1.32(s),4.72(h) 






44.4+48.5 
















Hard src 

Hard src 


P 
B 


14.3^13^^ 
10.515.6 


, ofi+1.94 
1-8^-0.46 


„ „7+0.67 
^-'5^-0.58 




2 12+0.80 
''-1^-0.42 


6 46+0-52 
°- °-0.46 


33, 

34, 


.1/28 

,5/28 


1 03+0-08 
1-O3-0.07 
1 02+0-08 
1-0''-0.07 


3 42+0-28 
''-*^-0.24 
2 07+0.16 
^-'"-0.14 


Total 


P 


2-47i°o:f7 


1 11+0.20 
-0.18 




0.0+0.02 


2 23+0.13 


20.0+0.06 


90, 


.0/28 


6.08+0.12 


11.81+0.30 


Total 


M 


1.33+0.27 


115.4+12.0 






64, 


.3/27 


5.85+0.63 


7.84+0.86 


Total 


PM 


0.90+0-4 


34+0-08 
" -0.07 


22.9+2.1 


0.23(bf) 


1 ^0+0.28 
1-50-0.27 


13.3+0.04 


21 


.9/26 


5.64+0.60 


6.90+0.76 



Notes to Table 1: 

^Source region (see paper 1). 

^Background accumulated from field for halo, soft sources, and total; background accumulated from diffuse halo and galaxy disk emission for disk, core, and 
hard point sources. For the PSPC spectral analysis we examine both "field" and "local" backgrounds. The held background is determined in source-free 
areas outside the galaxy that are selected to be larger than the source regions so as not to be background-noise limited, while the local background is 
determined in regions within the halo and galaxy disk that contain diffuse emission but that are free of compact sources. The latter is important within the 
central regions where the spatially-resolved spectra of the galaxy disk and nucleus can be contaminated by soft, diffuse emission that is seen in projection, 
"^Models are (P) power law, (M) MEKAL plasma, and (B) thermal bremsstrahlung. All models include absorption. 

^Absorbing column density in units of 10^^ cm~^. Where two values are listed, the first is for the soft component and the second is for the hard component. 

Temperature of thermal component in units of keV. Where two values are listed the first is for the soft component and the second is for the hard component. 
For all fits, the absorbing column was constrained to be greater than the Galactic value of 9 X 10"'"^ cm ^. 

■^Normalization of thermal component in units of lO^'^K, where K is 10~^'^/(4-7r D*^) J nenjjdV for the MEKAL model and 3.02 x 10~^^ /(47r D^) 
J nenidV for the bremsstrahlung model. D is the luminosity distance to the source in cm and ne, npj, and nj are the electron, hydrogen, and ion densities 
in cm ^ . 

^Abundance in units of Zq. 

Power-law photon index. 
*Power-law normalization in units of 10""* photons keV"''^ cm~^ s"''^ at 1 keV. 
^ divided by the number of degrees of freedom. 

''^Observed and Unabsorbed 0.1 — 2.0 keV fiuxes in units of 10~^^ ergs cm~^ s~^. Where two values are listed, the first is for the soft component and the 
second is for the hard component. 

Except for fluxes and normalizations, which have 90% confidence errors for one interesting parameter (x^ + 2.71), errors represent 90% confidence errors 

for 2 or 3 interesting parameters (x^ + 4.61 and + 6.25). 
"(f)" denotes a fixed parameter. 

"(bf)" denotes a best-fit value that is fixed to determine errors on the other parameters. 
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Table 2: Comparison of ASCA modeling for NGC 253* 



Model" 
(1) 



norm" N-^(w)^ 



(2) 



(3) 



1.9 
+0.9 
1.0 



kT(w)" 
(4) 



(5) 



OMSS 
(6) 



-i.U4 
'■""-0.70 
0.21±0.08 



m 



11.7+6-0 

„+3.l' 



r/kT^ 
(8) 



2.„7l0.22 

7 01+2-74 
- -1.67 

«-S^ll.69 



(9) 



(10) 



t(c)9 
(11) 



f(w)9 
(12) 



f(h)9 
(13) 



PMM 
PM 
BMM 
BM 



4.0 
1.5 



4.6"' 



R+1-6 
-1.5 

, q + 0.9 

^•^-0.9 



0.60±0.06 
0.64±0.07 
0.61±0.06 
0.65±0.08 



1.55 
10.7 
1.68 
9.47 



2.27"" 



7+1.24 
-0.79 

°-^^l°o:^^7 



5.0 



2.3 
+4.0 
-3.3 



2.04 
2.28 
1.59 
1.66 



515.4/456 
525.9/457 
626.0/456 
640.0/457 



1.29 1.00 

2.23 0.71 

1.23 1.07 

2.19 0.78 



Notes for Table 2: 

^Models are M: Mewe-Kaastra plasma (MEKAL), P: power law, B: bremsstrahlung. We assume 2j^p = 2qj.^ = 1.0 for 3c and z?[.^p = Z(^j.^=0.05 for two- 
component models, as found by P97. For the cool component kT = 0.2 keV and the abundance is 0.2 times solar (see joint fits). The Galactic column 
density is included for all fits. 
^Normalizations are given in units of 10~^. 
°The units on JVjj are 10^^ cm~^. 

^Temperature of the warm component in units of keV. 
^Abundance of 0,Mg,Si, and S in solar units. 

^Photon index (temperature in keV) of the hard (hot) component. 

^Fluxes arc 0.1-2.0 keV fluxes in units of 10 12 ergs cm ^ g 1 

*The spectral extraction regions used here are the same as those of P97. 



Ta ble 3: NGC 253 ASCA and PSPC Parameters from Joint Fits 

Soft Medium Hard 



Model" kT^ JVj^(w) fcT| Zhcn ^qmSS ^NF F/kT x^/i' 
O) (2) (3) (4) (5) (6) (7) (8) (9) (10) 



PMMl 


0.26+0.04 


0. 


18+0-11 
^^-0.05 


-«1^3;^^ 


0.70+0.10 


0.5(f) 


o-«li;^ 


0-^li;^ 


8 51+3-E>4 
®-^^-3.05 


1.96+0- 


18 
,20 


538.5/499 


PMM2 


0.26+0.04 


0. 


17+0.11 
-0.05 


6.77^3.03 


0.72+0.09 


0.5(f) 


o-srlJ:f7 


0.6(f) 


7 Ql +3.99 
'-'^ -3.26 


1.96^0; 


17 

19 


640.5/500 


BMM 


0.27+0.04 


0. 


17+0.10 

-0.05 


7 67+3-23 
'-""-4.62 


0.75 + 0.09 


0.5(f) 


96+2-34 
-0.58 


0.5(f) 


4 60+''-30 
*-''"-2.63 


7.8t^; 


7 



544.2/500 


PMM3 


0.27+0.03 


0. 


, 7+0.09 
-0.05 


74+0.16 
-0.14 


0.5(f) 


79 + l-''2 
"-'^-0.45 


0.5(f) 


=NH(m) 


1.9210; 


14 
16 


540.8/501 


PMM4 


0.27+0.03 


0. 


,7+0.09 
-0.05 


7.28t^-« 


73+0-11 


1.0(f) 


1 67+0-67 
-°'-0.75 


1.0(f) 


=NH(m) 


1.94i0. 


10 
11 


540.6/501 


(PMM)pc 


0.34+0.07 




Z{hard) 




72+0-30 


1.0(f) 


o,+3.91 
■^-^°-2.17 


l.O(t) 


*12 Q+4-3 


2.12+0- 


18 
,17 


640.5/500 



Notes to Table 3: 

'^Model components are P:power law, B: bremsstrahlung, M: MEKAL plasma, pc: partial covering. 
^Temperature of cool component in units of keV. 
'^Abundance of cool component in units of Zq. 

'^Column density of warm component in units of 10^''" cm~^. 
^Temperature of warm component in units of keV. 
^ Column density of hard component in units of 10'^"'^ cm""^. 
Errors for NH(w), NHi 
both columns are free. 

The Galactic column of 0.9x10^^ included in all fits. 

*For the partial covering model, the covering factor is fc = '^■^^^q^3 
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Table 4: NGC 253 Normalizations and fluxes from Joint Fits 



Model 


A(c)'' 


A(w)« 


A(h)« 






{{hf 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


PMMl 
PMM2 
BMM 
PMM3 


C qq + U.8b 

c 40+O.8O 

c: «1+0.75 
^•"J^-0.89 

^•"3°-0.19 


4.27+i;^^ 
3.64tl;f^ 

q 79+0.85 
•^^ ' ^-0.67 

4.UJ_Q 74 


■^■^^-0.29 

1 62+"-" 
r-yu_o.i9 


2.93 
2.96 
2.99 
3.00 


1.04 
0.95 
0.87 
0.90 


1.04 
1.09 
1.15 
1.09 


PMM4 


c C4+O.2I 
^•^^-0.20 




1.98±0.15 


2.98 


0.86 


1.16 


(PMM)pc 


0.34±0.10 


0.78t«;ll 


2.69±0.33 


1.60 


0.88 


3.17 



Notes to Table 4: 

"Normalizations for the cool, warm, and hard components are given in units of 10 
^Fluxes are 0.1-2.0 keV fluxes in units of 10 ergs cm ^ s ^. 



Table 5: PSPC Spectral Results for M82 



(1) 

Region 

(1) 


(2) 

M,F 

(2) 


(3) 

JVh 
[1020 cm-2] 
(3) 


(4) 
kT 
[keV] 
(4) 


(S) 

^kT 
[10-3] 
(5) 


(6) 

Z/Zq 

(6) 


(7) 

r 

(7) 


(8) 

Ar 

[10-3] 
(8) 


(9) 
(9) 


(10) 
fx 
[Obs.] 
(10) 


(11) 

fx 
[Unabs.] 

(11) 


(12) 

Lx 
[Unabs.] 
(12) 


Halo 
Halo 


P 
M 


'■"-11.0 
7 5+1.5 

'■'^-i.e 


0.45+0.06 


8.10+0.40 


0.03+0.02 


6.32+-", 


^- -0.17 


77.0/27 
34.9/26 


-«l"o:i^ 

1.68+0.07 


4.27+0.21 


6.42+0.26 


Halo 


MM 


4 4+2-0 

-0.7 


0-3^-0.05 


0.67+0.05 


1.0(f) 






14.3/25 


1.60+0.16 


2.26+0.19 


2.87+0.26 



Core 
Core 
Core 
Total 
Total 
Total 



P 
PMl 
PM2 

P 

M 
PM 



33.5^ 



^+3.5 
-2.6 
76+2-3 



31.8 
31.4 
16.4 



1.3 
12.1 
1.6 



+ 3.2 
3.3 
4.7 
1.4 
7+0.8 



r,-{-unb 
-2.82 



67+°-04 
-0.05 

54+0-10 
"■° -0.16 



97+0-08 
■ -0.17 

60+0-05 
"°"-0.03 



, Q2+0.13 

1 92+1-15 
l-''2_o.45 



■'^■"'-4.0 
2 ns+0-13 
2-08_o,i4 



1.0(f) 

1.0(f) 
1.0(f) 

0.02+0.01 
1.0(f) 



2 qa+0-19 
^-''*'-0.18 
„ ,„+0.26 
•^^-0.29 

1.50(f) 

o ,r + 0.2 



7+0.09 
-0.08 



^+0.2 
-0.1 
„ Qo+0.27 
2-98_o.26 

+0.33 
0.26 
H.4 
-1.1 



12.5"' 



5.37 



3 34+0-35 
"■"^-0.33 



502.1/27 
41.2/25 

46.3/28 
568.3/166 
540.2/165 
124.7/164 



10.810:3 
11.1+1.0 

8 56+1-46 
■ -0.91 

12 4+1-3 
-1.3 
12 3+3-7 
-0.5 
H.O 



16.1+1.2 



19.18+1.52 
7+6.7 
-3.6 



33. 7"^ 



12.7"' 



1.6 



Notes to Table 5: 

Models, normalizations and errors are the same as in Tables 3 and 4. 

in Q 1 OQ 1 

Fluxes are 0.1—2.0 keV fluxes in units of [10 ergs cm s ]. Luminosities are 0.1—2.0 keV luminosities in units of [10 ergs s ]. 
The core spectrum has local background subtracted. 
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Table 6: Comparison of ASCA modeling for M 82 



Model" 


norm'' 


JVH(m) = 


kTCw)" 


norm^ 


Z 


iVH(h) = 


r/kT-'' 


norm'' 


X^/dof 


f(c)9 


f(w)9 


f(h)9 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 



PRs 

PMs 

PRv 

PMv 

BMv 

*PRRs 

*PMMs 



1.31±0.07 
1.53±0.07 



—4.5 
1 1+0-6 
-0.5 
1+0.5 
-0.3 
,+0.3 
-0.4 
1+0.3 
-0.4 

7.60l0;« 

8.40±0.9 



11_ 
2.2^ 
2.1^ 



0.77±0.03 
0.60+0.03 
0.75+0.05 
0.61+0.04 
0.61+0.03 
0.77+0.03 
0.67+0.04 



1.22+0.16 
2.70+0.11 
29.5+1.4 
32.3+2.2 
29.4+2.4 
12.6+0.3 
16.7+0.4 



1(f) 
1(f) 

(a) 
(b) 
(c) 

1(f) 
1(f) 



r,+3.6 

7.7 
+3.9 
3.3 
+3.0 
2.7 
g 9+1.8 

10.6+^-^ 



25.5 
16.7 
11.6 



1.50+0.04 
1.55+0.04 

1.77+0.08 
1.71+0.10 

1.47+0.08 
1.54+0.09 



3.68+0.06 
3.91+0.10 
6.19+0.10 
5.42+0.09 
5.44+0.08 
3.60+0.10 
3.96+0.06 



1205.0/813 
1167.7/813 
903.2/805 
879.3/804 
875.7/804 
969.6/811 
952.4/811 



3.0+0.2 
3.3+0.2 



8.5+0.2 
4.0+0.2 
11.1+0.6 
9.7+0.8 
9.4+0.8 
5.5+0.2 
5.3+0.2 



3.0+0.4 
6.6+0.3 
1.0+0.7 
1.5+0.1 
1.7+0.1 
2.3+0.1 
2.2+0.1 



Notes for Tabic 6: 

^Models arc R: Raymond-Smith plasma, P: power law, B: bremsstrahlung. 

For the cool component, the temperature is fixed at 0.31 keV and TVj^ — 3.71 X lO'^*^ cm"^ 

(a) abundances are; N= 1 (f) ,0 = 0.0, Nc=0.16, Mg=0.26, Si=0.35, S = 0.6(), Ar,Ca,Ni = 0.014, Fe=0.034 (Si/Fe = 10.3). 

(b) abundances arc: N = l{f), = 0.19, Ne=().2(), Na=l(f), Mg=0.33, A1=0.0, Si = 0.45, S = 0.S4, Ar,Ca,Ni=0.0, Fe=0.07 (Si/Fe = 6.4). 

(c) abundances are; N = l(f), = 0.21, Ne=0.23, Na=l(f), Mg=0.34, A1=0.0, Si=0.49, S = 0.97, Ar,Ca,Ni=0.0, Fe=0.08 (Si/Fe = 6.1). 
Fluxes are averaged GIS fluxes for M 82. 

= 6 X 10'^" cm~-^ and kT{c) = 0.35 keV. 



Table 7: M82 ASCA and PSPC Parameters from Joint Fits 







Soft 






Medium 




Hard 




Model" 
(1) 


JV^(c) 
(2) 


kT\ 
(3) 


(4) 


JVj^(w) 
(5) 


(6) 


(7) 


iV^j(h) 
(8) 


F/kT 
(9) 


(10) 


PRs 
PRv 

PMs 


37+1.5 

4 6+11 
., O+0.6 

■^■**-o.i 








"■'°-0.03 
0.69+0.03 

0.64+0.02 


1.0(f) 

(a) 
1.0(f) 


^■l+53 

^■s^iSi? 

14. ±3. 9 


1 r,4 + ().U3 

■ -0.04 
1 44+0-()-'i 

■ -0.06 


1396.3/975 
1095.1/968 

1360.9/975 


PMv 


5 72+1-2 






„ -+0.6 
^■^-0.7 


0.60+0.03 


(b) 


2 4+"-9 
-0.8 


l-«in? 


1066.6/968 


PMMs 


5 3+1-" 


0.33+0.03 


1.0(f) 


0.68+0.05 


1.0(f) 


10.113.6 


, _,+0.09 
^■^3-0.08 


1110.9/972 


PMMv 


6 2+1-0 


0.48+0.05 


(d) 


4 5+1.6 
■°-1.4 


"■''''-0.12 


(d) 


6 l+l-'' 
"■ -1.8 


1.47+0.09 


1045.5/965 



Notes to Table 7: 

"Model components are P: power law, R: Raymond- Smith plasma, M: MEKAL plasma. The s and v subscripts refer to solar and variable abundances, 

respectively. 

''Temperature of cool component in units of keV. 

"^Abundance in units of .2q . 

Column density of warm component in units of 10 cm . Column density of soft component in units of 10 cm 

Temperature of warm component in units of keV. 
^ Column density of hard component in units of It)^"'" cm ^. 

Errors for NH(w), NH(h), V . and Z arc 90% confidence for 4 interesting parameters (Ax*^— 7.1). Errors for kT{w) are Ax —4.6 for cases where both columns 
are free. 

The Galactic column is 3.7x10*^'-' cm"'^. 

(a) abundances are: N=l(f), = 3.17, Ne^4.47, Mg^6.41, Si^l3.6, S^33.8, Ar,Ca,Ni^S.93, Fe^O.88. 

(b) abundances are: N=l(f), 0=4.34, Ne-5.63, Na-l(f), Mg-9.43, Si-15.5, S-29.5, Ar,Ca,Ni=4.88, Fe=1.98. 

(c) abundances are: N=l(f), = 1.94, Ne^4.70, Mg^3.97, Si^3.58, S^5.26, Ar,Ca,Ni^8.98, Fe=2.10. 

(d) abundances are: N=l(f), 0=0.85, Ne=0.01, Na=l(f), Mg=1.78, Si=2.05, S=2.41, Ar,Ca,Ni=3.96, Fe=0.51. 
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Table 8: M 82 Normalizations and fluxes from Joint Fits 



Model 


A(c)" 


A(w)° 


A(h)" 


f(c)'' 


f(w)'' 




PRs 




1.39±0.08 


4.98±0.05 




3.5 


8.9 


PRv 




1.11±0.05 


4.24±0.04 




5.6 


5.5 


PMs 




1.75±0.07 


5.36±0.05 




4.1 


8.2 


PMv 




1.02±0.04 


4.54±0.04 




6.2 


5.1 


PMMs 


2.35±0.10 


21.1±0.6 


4.92±0.07 


4.8 


5.9 


2.9 


PMMv 


3.78±0.15 


5.54±0.14 


4.54±0.06 


6.1 


2.9 


3.4 



Notes to Table 8: 

"Normalizations for the cool, warm, and hard components are given in units of 10 
''Fluxes are 0.1 — 2.0 keV observed fluxes in units of 10^^^ ergs cm~^ s~^. 



Table 9: Spectral Models for Spatially-Resolved Emission 



Galaxy 


Halo 


Disk 


Core 


Pt. 


sources 


Notes 


(1) 


(2) 


(3) 


(4) 




(5) 


(6) 


NGC253 


M+M 


M+M 


P+M 


P 


or M 


absorption in disk 


M82 


M+M or P+M 




P+M 


P 


or M 


absorption in disk and halo 



Notes to Table 9: 

Models are M: Mekal plasma, P: power law. 
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90 

Fig. 1. — MEKAL plasma model for a kT = 0.5 keV solar-abundance plasma (top curve). The strongest 
emission lines from each element are labeled with the exception of the Fe and Ne complexes, which are 
shown separately for clarity in the middle and bottom curves. The Y axis is arbitrary. 



23 




Fig. 2.— The Rosat HRI image of NGC 253. The field is 25' x 25' and the data are smoothed for an 
effective resolution of ICK'5 (of. paper 1). 
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Fig. 3. — PSPC image binned to have 15" pixels and consisting of photons that fall within the 0.25 keV 
band and 1.50 keV band as defined in paper 1. The circular and elliptical regions used to extract the 
spatially-resolved spectra of the galaxy halo, disk and point sources are drawn. 
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Fig. 4. — PSPC spatially-resolved spectra of the diffuse disk emission with compact sources removed. For 
(a) the field background is subtracted and for (b) the local background is subtracted. Crosses represent 
the southern disk and triangles represent the northern disk. 
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Fig. 5. — Model-fitting results for NGC253. Open squares (disk emission), diamonds (halo, soft point 
sources), stars (core) and solid triangles (hard point sources) are the spatially-resolved PSPC results for 
single-component model fits (either a power law or a MEKAL plasma), except for core(2), which results 
from a two-component power law plus MEKAL plasma model fit to the core spectrum. Open triangles 
are the ASCA result from the corresponding spectral component from P97. Circles indicate that a fit is 
statistically unacceptable as defined by xt greater than 1.3. (a) A'^h vs. F for a power-law model. The 
point labeled P97 represents the corresponding hard power-law component from P97; core(2) represents 
the hard core component from the PSPC (model PM3, Table 1). (b) Ah vs. kT for a single MEKAL model 
with free abundances. The point labeled P97 represents the corresponding thermal component from P97; 
core(2) represents the soft core component from the PSPC (model PM3, Table 1). For comparison we also 
plot the single-component model result for the integral PSPC spectrum (labeled "Total"). Spectra of the 
halo and soft compact sources have field background subtracted; those for the disk have halo background 
subtracted; those for the core have halo plus diffuse disk emission subtracted. 
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Fig. 6. — Our best fitting ASCA model for NGC 253 which consists of a two-temperature MEKAL plasma 
with temperatures of kT = 0.2 keV and kT = 0.6 kcV and a power law with r=2.0. The abundances are 
fixed at 0.2 Zq for the cool component and 1.0 Zq for all but O, Mg, Si, and S for the warm component. 
The derived absorption is 4 x 10^^ and 84 x lO'^^ cm^'^ for the warm plasma and power-law components, 
respectively; see Table 2 (model PMM). 
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Fig. 7. — (left) ASCA and Rosat data for NGC 253 and the best-fitting two-temperature MEKAL plus 
power-law model folded through the instrumental response (Table 3, x^=540.5). The ratio of the data 
to this model is plotted in the bottom panel, (right) The model (Table 3 model PMM2). 




Fig. 8. — (left) ASCA and Rosat data for M 82 and the best-fitting two-temperature MEKAL plus power- 
law model folded through the instrumental response (see Table 7, model PMMy). The ratio of the data 
to this model is shown in the bottom panel. A Gaussian to model the Fe K emission at ~ 6.5 keV is not 
included in the model; the line is visible in the ratio plot. 
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